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Abstract 

The  transport  phenomena  of  mass  and  electrical  charges  play  an  important  role  in  the  proton  exchange  membrane  utilised  in  fuel  cell 
applications.  The  calculation  of  the  membrane  ion  and  water  transfer  is  usually  determined  in  two  ways:  one  is  a  semi-empirical  mode  that 
can  be  experimentally  set  for  any  application  and  cannot  be  used  in  the  prediction  of  the  characteristics  of  new  cells,  and  the  other  is  a  very 
complex  mathematical  approach  that  needs  very  long  calculation  times.  In  fact,  the  membrane  mass  transfer  model  can  cause  problems 
in  the  management  of  PEM  stack  models,  by  increasing  calculation  times  and,  above  all,  convergence  problems.  This  work  attempts  to 
overcome  this  bottleneck  using  a  novel  approach.  The  validation  of  the  new  approach  has  been  made  by  applying  the  model  results  to  a 
PEM  fuel  cell  model  and  by  comparing  the  overall  results  to  those  in  the  literature.  The  resulting  calculation  time  warranted  the  application 
of  this  model  to  a  very  complex  stack  simulation,  resulting  in  good  fit  and  reasonable  computing  times. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  fuel  cell  the  direct  conversion  of  the  chemical  energy 
of  the  reactants  to  electrical  energy  is  achieved  with  high  ef¬ 
ficiency  and  high  environmental  compatibility  [1].  Among 
the  different  types  of  fuel  cells,  proton  exchange  membrane 
fuel  cells  (PEMFC)  are  now  considered  the  best  solution  for 
automotive  applications  [1-7].  This  is  due  to  the  flexibility 
of  their  use  (starting  is  very  simple,  the  device  is  of  small 
dimensions,  efficiency  is  very  good)  and  their  low  work¬ 
ing  temperature,  which  is  important  for  the  intrinsic  safety 
of  the  vehicle.  For  these  applications,  the  modelisation  of 
the  stack  is  needed  to  proceed  to  the  definition  of  the  over¬ 
all  vehicle  design.  In  particular  [5]  the  use  of  a  dynamic 
model  can  be  considered  the  basis  for  the  successful  assem¬ 
bling  of  all  the  power  devices  (stacks,  feeder,  batteries,  util- 
isators).  The  dynamic  response  can  be  pursued  only  if  the 
stack  model  does  not  need  high  calculation  times  and  in  this 
case  some  simplifications  may  be  necessary  [5]  in  the  cal¬ 
culation  of  membrane  conductivity,  an  intrinsic  property  of 
PEMFCs. 
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These  cells  consist  of  a  perfluosulfonic  membrane  elec¬ 
trolyte  coupled  with  two  electrodes  (anode  and  cathode) 
made  of  a  substrate  of  carbon-PTFE  coated  with  platinum 
as  electrocatalyst,  as  schematised  in  Fig.  1. 

Air  is  fed  to  the  cathodic  compartment  while  hydrogen  is 
fed  to  the  anodic  one  and  the  electrolyte  performs  the  func¬ 
tions  of  transferring  H+  from  the  anode  to  the  cathode  and 
separating  gaseous  reactant.  In  particular,  the  H+  transfer 
can  be  realised  only  when  the  membrane  is  strongly  hy¬ 
drated  [5]. 

Electrochemical  reactions  involved  are: 

H2  -+  2H+  +  2e_  (anode) 

2H+  +  \02  +  2e_  -+  H2O  (cathode) 

H2  +  +2  -*  H20  (overall) 

Stacks  are  obtained  by  connecting  a  number  of  cells  in  series 
with  gas  distributors  (bipolar  plates)  that  feed  the  fuel  to  the 
anodic  side  and  the  oxidant  to  the  cathodic  one. 

As  mentioned  above,  hydration  level  conditions  the  pro¬ 
ton  transfer  (or  proton  conductivity)  in  the  membrane.  In 
fact,  the  sulfonic  groups  easily  dissociate  into  SO3-  and 
H+  in  these  membranes  in  the  presence  of  water.  Under  this 
condition,  the  proton  can  be  considered  as  a  mobile  charge 
that  encounters  a  low  resistance  when  moving  across  a  po¬ 
tential  gradient.  To  achieve  the  membrane  humidification, 
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Nomenclature 

a 

air  relative  humidity,  dimensionless 

Ci 

concentration  of  (molm-3) 

Ci  s 

dry  ionic  concentration  of  “i”  (molm-3) 

Di 

diffusivity  of  in  the  membrane 

(m2  s-1) 

Dx 

diffusivity  in  the  hydrated  membrane 
defined  by  Eq.  (2)  (m2  s-1) 

e 

humidification  parameter,  dimensionless 

F 

Faraday  constant  (Cmol-1) 

I 

electrical  current  density  (Am-2) 

j 

adimensional  current,  defined  by 

Eq.  (21) 

J 

total  flux  of  ions  (molm-2  s-1) 

Ji 

flux  of  (molm-2  s-1) 

k 

hydraulic  permeability  of  the  membrane 
(m2) 

h 

Boltzmann  constant  (Jmol-1  K-1) 

m 

membrane  thickness  (m) 

M 

equivalent  weight  of  membrane,  i.e.  weight 
of  membrane/mol  of  SO3-  groups 
(g  mol-1) 

n 

adimensional  total  water  flux,  defined 
by  Eq.  (21) 

nc 

number  of  electrons  transferred  in  the 
reaction,  dimensionless 

P 

total  pressure  (Pa) 

R 

gas  constant  (Jmol-1  K-1) 

s 

membrane  thickness  in  dry  conditions  (m) 

T 

temperature  (K) 

V 

potential  (V) 

Vi 

molar  volume  of  (m3  mol-1) 

X 

geometrical  coordinate  along  the 
membrane  (m) 

geometrical  coordinate  along  the  dry 
membrane  (m) 

y 

adimensional  geometrical  coordinate 
defined  by  Eq.  (17) 

yi 

molar  fraction  of  “2” 

Zi 

electric  charge  of  “2” 

Greek  letters 

a 

parameter  defined  by  Eq.  (30), 
dimensionless 

p,  p',  p" 

parameters  defined  by  Eq.  (22), 
dimensionless 

y 

yi,  y 2, 

parameter  defined  by  Eq.  (21), 
dimensionless 

y  3,  Y4 

coefficients  in  Eq.  (3) 

r'v  Y2 

coefficients  in  Eq.  (32) 

8,  8' 

parameters  defined  by  Eq.  (22), 
dimensionless 

K 

electro-osmotic  drag  coefficient, 
dimensionless 

T) 

adimensional  membrane  tension 

0 

diameter  of  water  molecule  (m) 

X 

membrane  hydration,  i.e.  the  ratio 
H2O/SO3-  in  membrane,  dimensionless 

TC 

adimensional  pressure 

P 

membrane  density  (kgm-3) 

O 

<Pl,  <P2,  ^3, 

resistivity  m-1) 

<P4,  <P5,  <P6 

coefficients  in  Eq.  (2) 

X 

ratio  between  the  net  water  flux  through 
the  membrane  and  the  water  produced 
by  the  electrochemical  reaction, 
dimensionless 

Xip->  Xiv 

f  l,  f2. 

non-ideality  coefficients,  dimensionless 

^  3,  if  a 

coefficients  in  Eq.  (27),  dimensionless 

Q 

electrical  resistance  (£2  m2) 

CO 

adimensional  resistance 

Subscripts 

A 

referring  to  the  SO3-  groups  in  the 
Nation®  membrane 

H 

hydrogen 

old 

old  value 

W 

water 

Wt 

total  water 

0 

anodic  side 

1 

cathodic  side 

Superscript 

* 

reference  condition 
(i.e.  complete  hydration) 

reactant  gases  are  usually  humidified  by  being  bubbled 
through  high-temperature  water  columns  before  being  fed 
into  the  cell.  Flooding  and  dehydration  of  the  membrane  are 
phenomena  to  be  avoided:  when  the  first  situation  occurs 
the  pores  of  the  electrodes  and  the  bipolar  plates  can  be 
obstructed  by  water  and  will  not  be  accessible  to  the  gases; 
in  the  second  case  the  membrane  conductivity  decreases  to 
zero  and  also  the  damages  of  the  membrane  can  occur. 

The  literature  on  PEMFC  modelling  includes  studies  on 
the  local  electrochemical  kinetics  of  the  membrane-electrode 
assembly  [8-12]  based  on  the  analysis  of  mass  and  ion 
transfer.  However  the  application  of  these  approaches  [13] 
involves  long  calculation  times  for  only  a  point  in  the  I-V 
cell  characteristic  curve. 

In  order  to  set  up  a  tool  for  the  dynamic  prediction  of 
cell  behaviour  and  avoid  the  problems  involved  in  other 
approaches,  a  novel  approach  to  the  problem  is  presented  in 
this  paper. 
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2.  Mass  and  charge  fluxes  in  conductive  polymeric 
membrane 

In  a  PEM  fuel  cell,  proton  conduction  through  the  mem¬ 
brane  is  substantially  due  to  the  transport  of  protons  in  the 
solid  electrolyte.  The  most  used  polymer  for  this  applica¬ 
tion  is  Nation®,  which  has  the  chemical  structure  shown  in 
Fig.  2. 

The  most  important  properties  of  Nation®  for  ionic  con¬ 
duction  in  PEMs  are: 


they  are  formed,  to  that  cathodic  one  where  they  will  meet 
the  oxygen  and  react). 

•  Water  diffusion,  that  moves  from  the  cathode  side,  where 
the  water  is  produced,  to  the  anodic  one  where  a  small 
concentration  occurs.  Hydrogen  ions  are  involved  in  a 
counter-diffusion  mechanism  (from  the  anode  to  the  cath¬ 
ode). 

•  Pressure-driving  transfer,  which  happens  when  a  pressure 
difference  is  created  between  the  anodic  and  the  cathodic 
zones. 


•  the  presence  of  negatively  charged  sites; 

•  the  ability  to  be  hydrated; 

•  the  ability  to  transport  protons  between  the  charged  sites 
with  a  “jump”  mechanism,  thanks  to  the  presence  of  water. 

A  great  number  of  works  [9,12-17]  have  experimentally 
underlined  that  the  membrane  conducts  only  when  hydrated. 
The  hydration  originates  in  the  humidity  of  the  reagents  and 
in  water  produced  by  the  cathodic  reaction.  Indeed,  some 
phenomena  establish  a  water  concentration  gradient  in  the 
membrane. 

The  ion  transfer  phenomenon  can  substantially  be  de¬ 
scribed  together  with  other  phenomena,  such  as: 

•  Electro-osmotic  water  transport  from  the  anodic  side  to 
the  cathodic  one  (due  to  water  “dragging”  by  the  hydrated 
protons  that  are  transferred  from  the  anodic  side,  where 


mobile  proton 


Fig.  2.  The  chemical  structure  of  the  Nation®  polymer. 


As  can  be  seen  in  Fig.  1  the  ionic  conduction  of  the  mem¬ 
brane  can  be  described  by  the  attributes  of  the  charged 
species  and  of  the  water.  Table  1  underlines  the  species  name, 
its  charge  and  mobility  expressed  in  terms  of  diffusivity. 

According  to  the  semi-empirical  approach  of  Springer 
et  al.  [17],  extended  by  Costamagna  [13],  the  membrane 
overall  mass  balance  assumes  the  following  expression: 

I  I  p  dX  p  3n0  P\  —  P0 

X -  =  fr - —  Dk - kX  —  Dx - - - -  (1) 

nQF  nQF  M  ox  M  kbT  m 

The  term  on  the  left-hand  side  of  Eq.  (1)  represents  the 
total  water  flux  through  the  membrane;  the  terms  on  the 
right-hand  side  are  respectively  the  electro-osmotic  drag,  the 
water  diffusion  and  the  pressure  gradient  driven  flow  for  the 
water  in  the  pores  of  the  membrane.  The  adimensional  ratio 
between  the  membrane  water  concentration  in  the  membrane 
and  the  charged  sites  (SCE-),  here  identified  as  A.,  is  usually 
called  “membrane  hydration”. 


Table  1 


Electro-active  chemical  species  present  in  the  membrane 


Species 

Name 

Charge  (z) 

Diffusivity 

so3- 

A 

-1 

0 

H+  +  £H20 

H 

1 

#0 

h2o 

W 

0 

#0 
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Fig.  3.  Measured  intradiffusion  coefficient,  Dx,  adapted  from  [17]. 


In  equation  [1],  the  diffusivity  has  been  found  [17]  to  be 
strongly  dependent  on  membrane  hydration: 


Dx  =  (jp\  +  (P2^  +  (P3^2  +  (P4^3)  exp 


<P5 


1 

<P6 


(2) 


Eq.  (2)  represents  an  “equivalent”  diffusivity,  measured 
when  the  membrane  is  hydrated  and,  consequently,  swollen 
by  the  presence  of  the  water  in  its  molecular  structure.  It 
does  not  accurately  represent  the  diffusivity  for  low  mem¬ 
brane  hydration  values  (i.e.  X  <  3),  as  can  be  seen  in  Fig.  3. 

The  membrane  conductivity  is  evaluated  in  [17]  as  a  func¬ 
tion  of  the  hydration  along  the  membrane: 


o  =  (M 


Y2)  exp 


V3 


So,  an  electrical  resistance  is  associated  to  the  membrane 
that  causes  a  membrane  voltage  loss,  which  assumes  the 
general  expression  of: 

r  dx 

y]  =  m,  q  =  /  —  (4) 

Jo  O 

Notwithstanding  this  simple  explanation,  this  model  as¬ 
sumes  without  explicit  mention  a  number  of  simplifications. 
In  addition  the  integration  of  Eq.  (1),  where  the  diffusiv¬ 
ity  is  expressed  by  Eq.  (2)  as  a  function  of  the  membrane 
hydration,  can  cause  a  long  calculation  time  and  some 
convergence  problems. 

In  the  following  pages  the  physical  explanation  of  Eq.  (1) 
will  be  focused  on  in  terms  of  the  flux  of  the  species,  a  robust 
convergence  criterion  will  be  proposed  for  the  integration 
and  an  approximate  analytical  solution  that  requires  short 
computing  times  will  be  pursued  by  introducing  a  number 
of  weighted  approximations. 

The  flux  of  the  species,  when  concentration,  pressure  and 
electric  potential  gradients  occur,  can  be  written  following 
the  equation  (see  Fig.  4): 


li  =  -A  (V  •  Ci  +  x.pff 


This  is  substantially  the  Nernst-Planck  equation,  which 
is  assumed  to  hold  along  the  membrane  thickness: 
non-linearity  effects  due  to  the  membrane  structure  and 
variable  local  hydration  can  be  taken  into  account  by  means 
of  non-constant  diffusivity  coefficients  D/(A)  and  correction 
factors  Xip.  Xiv- 

In  the  membrane  the  general  law  of  electroneutrality  can 
be  written  as: 

E  Oz,  =  0  (6) 

i 

By  applying  Eq.  (5)  the  total  mass  flux  and  total  current  is 
written  as: 

j  =  J2ziJ-i  (?) 

i 

1  =  JF  (8) 

Considering  only  the  gradients  on  the  x-axis,  the  flows  of 
the  three  species  are  written  by  applying  Eq.  (5) 


7a  =  0 


Anode  Membrane  Cathode 


Elcctro-osmotic  transport 


Diffusion 


Pressure-dri  venflow 


Fig.  4.  Transport  phenomena  in  the  Nation®  membrane. 
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Jr  =  —Dr 


dCH  CuVudP  F  dV 

+  XHp  ^  - - f  ZrCrXRv 


dx 


RT  dx 


RT  dx 


Jw  =  — 


dCw  Cw  Vw  dP 

£>w  I  — +  XwP 


dx 


RT  dx 


(10) 


(11) 


Avt  =  £X7  — 


Dy^C  As  Dw 


dCw  = 


CAs 


dCu  =  - 


1  +  e 

Cas 


1  - 


X  de 


1  ft-  e  dX 


dX, 


de 


(1  +  e)2  dX 


dX 


(19) 


Incorporating  Eqs.  (10),  (11),  (15)-(20)  in  Eqs.  (13)  and  (14) 
we  obtain  the  following  system: 


1 


J  = 


Dr  DwC as 


D* 


La  +  <o2 

1  de  dX 


1  - 


X  de 


dX  Vw 
+  ~r~xr  XWp 


X  d7T 


1  +  e  dX  )  dy 
XHp  Vh  d7T  XHv 


V* 

Kw 


(1  +  e )2  dy  _ 


dr] 


s  |_(1  +  e)3  dX  dy  (1  +  e)2  dy  (l+^)2dy_ 


(20) 


The  sulfonic-ion  flux  is  zero  because  this  ion  is  a  part  of  the 
molecular  structure  of  the  Nation®. 

Using  Eqs.  (6)  and  (8)  we  have: 


Where  D ^  is  the  diffusivity  of  water  in  a  reference  condi¬ 
tion,  namely  dry  conditions  or,  better,  complete  hydration. 
Simplifying  system  (20)  with  the  relationships: 


Caza  +  Crzr  =  0, 

Ch  =  Ca 

(12) 

JwtS 

Js 

d7T 

(21) 

J  =  zrJr  =  Jr, 

I  =  JUF 

(13) 

1^1  - 

d*wCa7 

]  ~  D^Cas  ’ 

}  dy 

According  to  [17],  each  H  ion  drags  a  number  of  water 
molecules  which  is  approximately  proportional  to  hydra¬ 
tion.  So,  the  water  total  flux,  obtained  coupling  Eqs.  (10) 
and  (11),  is: 

Avt  =  £7rX  +  Jw  (14) 


Eqs.  (9),  (11)  and  (14)  underline  the  fact  that  the  transport 
phenomena  in  the  membrane  are  sufficiently  described  by 
the  total  water  and  H- species  fluxes. 

In  such  equations  electrical  and  pressure  forces  can  be 
suitably  expresses  in  adimensional  terms 


VF 

RT 

p_F w 

RT 


(15) 

(16) 


For  what  concerns  the  geometrical  coordinate  x,  one  should 
remember  that  the  membrane  thickness  m  increases  when 
humidification  occurs  and  it  is  possible  to  write  a  general 
swelling  law  dependent  on  the  hydration  as  dx/dxs  =  1  +  e, 
where  e  is  a  function  of  X.  So,  it  is  defined  a  membrane 
coordinate  system  (x)  and  a  reference  one  (xs).  Then,  the 
latter  can  be  referred  to  the  dry  thickness,  so  defining  a 
dimensionless  coordinate  y: 


dx 

s(  1  +  e) 


(17) 


By  applying  Eq.  (17),  the  concentration  of  the  “W”  and 
“H”  components  of  the  membrane  and  their  dependence  on 
hydration  X  can  be  written  in  terms  of  swelling  and  dry  ionic 
concentration  Cas* 


CAs 

1  +  e 


(18) 


5  = 


 Dw 

1 

X 

de " 

~  n* 
uw 

.(l+e)2 

(l+e)3  dk_ 

Dw  Vw 

1 

ft'  - 

 Du 

AVyp£>^  (1  +e)2’ 

P 

D* 

D\\  VH 

1 

ft" 

AnpD^V^(  l+e)2’ 

P 

—  XHv 

S'  = 


1  de 

(1  +e)3dX’ 
Dr  1 
£>w  d  +  d2 


(22) 


where  Eq.  (21)  define  terms  which  are  constant  when  y  and 
X  vary,  while  Eq.  (22)  define  terms  which  are  functions  of  X. 

Inserting  Eqs.  (21)  and  (22)  into  Eq.  (20)  and  rearranging, 
the  system  is  written  as: 


dX 

dy 


ySX  —  n), 


dX 

dy 


4 1  J-yS'  + 


a 


d?A 

dy) 

(23) 


The  first  equation  of  system  (23)  is  equivalent  to  Eq.  (1)  and 
their  terms  can  be  explained  in  the  same  way  as  those  in 
Eq.  (1).  So,  p  is  essentially  an  adimensional  diffusivity  co¬ 
efficient,  where  the  effects  of  the  term  in  brackets  of  the  first 
of  Eq.  (22)  is  expected  to  be  much  less  important  than  the 
dependence  of  Dw  on  X.  In  other  words,  the  first  of  Eq.  (22) 
puts  into  evidence  to  what  extent  the  effective  membrane 
diffusivity  depends  also  on  swelling;  empirical  expression 
such  as  Eq.  (2)  can  obviously  show  only  global  dependence 
on  hydration. 

It  is  also  worth  noting  that  the  dependence  on  X  of  pa¬ 
rameters  such  as  (3  and  /3'  is  also  linked  to  the  expression 
chosen  for  the  diffusive  fluxes,  that  is  Eqs.  (5),  (10),  (11)  or 
in  terms  of  molar  fractions,  such  as  A  =  —DiCj  dy;/dx. 

On  the  other  hand,  for  what  concerns  the  second  of 
Eq.  (23),  it  is  apparently  rather  different  from  the  empirical 
form  (3).  Obviously  these  differences  cannot  be  dramatic  in 
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practice,  but  they  underline  the  necessity  of  a  coherent  for¬ 
mulation  of  water  and  proton  transport  phenomena  in  front 
of  a  more  empirical  description.  A  simplified  and  direct 
path  is  to  consider  that  the  largely  most  important  force  in 
determining  ion  flux  is  the  electric  one:  in  fact,  of  the  three 
terms  in  the  right  hand  of  Eq.  (10)  the  third  is  normally 
one  order  of  magnitude  higher  than  the  others.  So  it  can  be 
assumed  as: 


=  8'  =  0, 


(24) 


and  ft'  corresponds  to  the  membrane  conductivity  of  Eq.  (3). 

Eq.  (23),  reorganised  in  their  full  extension,  can  be  written 
as 
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ySX  —  n 


PU  -  &'y) 

-  n) 


dX 


(25) 


Usually,  as  the  membrane  is  not  uniformly  nor  completely 
damped,  the  diffusivity  is  function  of  A.  Experimental  evi¬ 
dence  can  be  organised  in  terms  of  effective  diffusivity  co¬ 
efficients,  as  in  Eq.  (2).  So,  the  problem  of  calculating  the 
integrals  in  Eq.  (26)  in  a  number  of  real  cases  can  be  over¬ 
come  by  considering  the  equivalence  of  their  coefficients 
with  that  of  the  Springer  model  and  adopting  some  approx¬ 
imations.  Any  improvement  in  knowledge  on  the  behaviour 
of  coefficients  as  functions  of  humidity  and  temperature  can 
be  easily  incorporated  in  Eq.  (26).  On  the  other  hand,  for  an 
approximate  calculation  the  simplifications  adopted  by  De 
Francesco  and  Arato  [5]  can  be  assumed.  In  fact,  consider¬ 
ing  that  the  diffusion  coefficient  of  Eq.  (2)  has  been  eval¬ 
uated  by  fitting  experimental  values,  it  is  often  possible  to 
adopt  a  suitable  mean  value. 

For  instance,  if  the  first  of  Eq.  (22),  which  defines  the 
dimensionless  diffusivity  of  water,  is  considered  and  a  linear 
dependence  e  =  (const.  A)  <<C  1  of  swelling  on  hydration  is 
assumed  as  first  approximation,  we  have 


and  then  integrate  along  the  membrane  thickness,  from  y  = 
0  to  1: 
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The  evaluation  of  the  membrane  voltage  loss  can  proceed  by 
calculating  a  first  attempted  value  of  “ft”,  here  called  nw,  by 
considering  the  differential  written  in  Eq.  (23)  as  replaceable 
with  the  finite  difference  between  the  extreme  values  of  the 
hydration  of  the  membrane,  which  are  a  function  of  the 
external  relative  humidity: 


y  =  0;  A  =  A0;  rj  =  rjo,  y  =  1;  A  =  Ai;  rj  =  rji, 

X  =  \j/i  +  faa  +  ip^a2  +  \lf4a 3  (27) 


and  a  constant  mean  value  for  /3  can  be  assumed  within  the 
limits  in  which  a  constant  mean  value  for  D\  is  acceptable 
[5].  Similarly  constant  mean  values  can  be  used  for  the  co¬ 
efficients  8,  f . 

Thanks  to  this  approximation,  the  first  of  Eq.  (25)  can  be 
integrated  with  a  straightforward  analytical  path,  obtaining 


n  =  a 


Aq  + 


Ao  —  Ai 

exp(of  /  p)  —  1 


oi  =  U  ~  8y 


(30) 


This  approach  allows  the  determination  of  the  water  flux 
through  the  membrane  with  a  very  fast  calculation  and  with 
results  comparable  to  those  obtained  without  any  simplifi¬ 
cation  [5].  The  integration  gives  also,  if  of  interest,  the  hu¬ 
midity  field  along  the  membrane  thickness: 


Then,  by  adopting  a  Newton  method  for  convergence,  the 
next  value  of  “ft”  is: 


1  -  ft  (/V(£A  -  ySX  -  «om))  dx 

ft  —  ft  old  H - ~x~ - - - 

A0  {P/iKjX  -  ySX  -  n old)2)  dA 


(28) 


Proceeding  in  this  way  the  calculation  can  converge  without 
problems.  When  the  water  flux  n  is  known  the  voltage  losses 
can  be  easily  obtained  by  the  second  definite  integration  in 
Eq.  (26). 

It  is  possible  to  analyse  two  limit  cases:  if  the  membrane 
is  uniformly  damped  the  coefficients  can  be  considered  con¬ 
stant  (at  uniform  temperature);  in  the  case  of  the  membrane 
being  completely  dry,  the  flow  of  water  tends  to  zero,  and 
it  can  be  concluded  that  membrane  resistivity  is  very  high, 
according  to  the  typical  behaviour  of  Nation®  membranes 
and  their  total  inability  to  conduct  current  (or  protons)  when 
not  humidified. 


A  —  Aq  +  (Ai 


exp  (ay/P)  -  1 
exp  (a/ /3)  —  1 


(31) 


The  membrane  resistance  can  be  then  obtained  from  the 
integration  of  the  second  of  Eq.  (25);  for  an  approximate 
evaluation,  instead  of  Eq.  (25),  reference  can  be  made  to  the 
position  (24)  and,  moreover,  the  membrane  conductivity  can 
be  expressed  as  a  linear  function  of  A  in  a  similar  way  as 
Eq.  (3) 

P"  =  -  )4  (32) 


Remembering  Eq.  (31)  and  integrating  the  membrane  resis¬ 
tance  will  be 


m-11  a  +  /31n((yiAi  —  y2)/(Ki^o  —  X2)) 

CO  =  -  =  - 

j  yin  -  YiOL 


(33) 


where  a  and  n  are  given  by  Eq.  (30). 
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3.  Fuel  cell  local  model 

A  local  kinetic  model  of  membrane  and  adjacent  gas 
phases  has  been  elaborated  by  using  the  approach  presented 
above.  Starting  from  the  Nernst  potential,  the  cell  voltage  is 
expressed  by  applying  the  usual  definitions  of  the  voltage 
losses  as  reported  in  the  literature  [5]  and  by  applying  the 
model  explained  here  for  the  membrane  voltage  losses.  In 
the  model  the  water  transport  effects  outside  the  membrane 
are  neglected. 

As  a  first  application,  assuming  that  a  local  model  is  suffi¬ 
cient  to  describe  a  little  PEM  laboratory  single  cell,  the  char¬ 
acteristic  curves  of  such  a  cell  are  evaluated  for  average  cell 
temperature  T  =  343  K  and  gas  humidification  a  =  0.99. 
The  results  are  compared  with  experimental  data  obtained 
under  the  same  conditions  [5,12]  in  Fig.  5. 

As  usual,  the  current-voltage  curves  can  be  schematised 
in  three  different  regions:  (i)  semi-exponential,  where  the 
cell  works  in/with  activation  losses  at  the  cathode  at  very  low 
current  density;  (ii)  linear  behaviour,  determined  substan¬ 
tially  by  ohmic  losses;  (iii)  voltage  drop  at  high  current  den¬ 
sity,  due  to  mass  transport  limitations  in  the  membrane.  The 
model  shows  a  good  agreement  that  provides  a  validation  of 
the  membrane  model  and  of  the  deduced  kinetic  model:  in 
fact,  by  comparing  these  results  with  those  reported  in  the 
literature  for  the  non-analytical  model,  the  main  numeric 
differences  are  less  than  5%,  a  value  that  can  be  compared 
with  the  experimental  errors  and  does  not  influence  the 
good  behaviour  of  this  approach.  In  addition,  the  calculation 
times  are  dramatically  reduced  from  some  minutes  to  a  few 
seconds.  So,  this  approach  provides  a  new  opportunity  to 
build  dynamic  models  of  stacks  containing  the  calculation 
time. 

In  addition,  the  possibility  of  having  good  results  without 
spending  excessive  time  in  calculations  will  allow  the  imple¬ 


mentation  of  this  local  model  in  direct  methanol  fuel  cells 
(DMFC),  by  simply  considering  the  following  phenomena: 

•  the  crossover  of  the  methanol  in  the  membrane  [18]; 

•  the  presence  of  methanol  in  the  calculation  of  activation 

and  concentration  voltage  losses  [18]. 

4.  Conclusions 

At  the  moment  PEM  fuel  cells  offer  a  very  promising 
technology  for  distributed  and  automotive  power  generation. 
High  efficiency  and  reduced  environmental  impact  are  the 
most  important  features  of  this  technology. 

Model  analysis  is  very  important  for  the  design  and  sim¬ 
ulation  of  a  power  generator  and  the  approaches  utilised  are 
still  affected  by  mathematical  complexity  that  limits  their 
use  in  dynamic  simulations. 

The  possibility  of  overcoming  the  time-analysis  problem 
has  been  studied  in  this  work.  In  fact  it  has  been  shown  that 
one  of  the  most  difficult  aspects  of  PEM  fuel  cell  modeli- 
sation,  i.e.  membrane  conductivity,  can  be  solved  in  an  ana¬ 
lytical  way  that  will  allow  us  to  simulate  the  cell  kinetic  in 
a  very  short  calculation  time,  obtaining  a  good  agreement 
with  experimental  data. 

The  kinetic  model  presented  here  can  be  used  as  a  basis  for 
the  future  development  of  fuel  cell  mathematical  approaches 
and,  above  all,  in  the  prediction  of  dynamic  behaviour. 

The  results  obtained  demonstrated  that  by  considering  a 
constant  diffusion  coefficient  for  the  hydrogen  ions  the  be¬ 
haviour  of  the  cell  is  correctly  represented.  This  result  will 
allow  for  further  development:  the  implementation  of  the 
swelling  experimental  data  and  their  application  to  a  wide 
range  of  uses  and  direct  methanol  fuel  cells  (DMFC)  by 
considering  some  particular  phenomena  that  are  influenced 
by  the  presence  of  methanol  on  the  anodic  site. 


current  density  [A  m'2] 


Fig.  5.  Comparison  between  model  results  and  literature  experimental  data. 
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